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Рассмотрены двумерные дискретные ортогональные преобразования, которые позволяют значительно сократить время сжатия 
кадра видеопоследовательности и могут быть реализованы в виде аппаратного кодека. 


Цифровые изображения занимают все большую 
часть информационного мира. Развитие сети Ин- 
тернет, систем видеонаблюдения, прогресс в техно- 
логии производства цифровых камер, сканеров и 
принтеров привели к широкому использованию ци- 
фровых изображений. Отсюда постоянный интерес 
к улучшению алгоритмов сжатия данных, предста- 
вляющих изображения. Сложность алгоритмов, ис- 
пользуемых для компрессии изображений, неуклон- 
но растет - сказанное касается не только объема вы- 
числений, но и идейных основ построения алгорит- 
мов. В настоящее время существуют следующие се- 
мейства алгоритмов сжатия изображений [1—3]. 

• Алгоритмы сжатия изображений без потерь. 

• Вейвлет-сжатие (включая стандарт ^РЕС2000, 
основанный на дискретном вейвлет-преобразо- 
вании ОІѴТ). 

• Фрактальное сжатие. 

• Алгоритм сжатия .ІРЕО и его модификации. 

Сжатие изображений из чисто теоретической за- 
дачи стало повседневной практикой, что требует 
при его осуществлении постоянного внимания к 
возможностям реальной аппаратуры и экономиче- 
ским характеристикам, например, для систем, обес- 
печивающих хранение, воспроизведение и обработ- 
ку видеоинформации в мобильных сетях класса 30, 
в системах видеонаблюдения, связанных с компью- 
тером по низкоскоростной линии связи (см. науч- 
но-техническую программу «Исследования и разра- 
ботки по приоритетным направлениям развития 
науки и техники» на 2005-2006 гг. по приоритетно- 
му направлению «Информационно-телекоммуни- 
кационные системы» ИТ-13.4/002). Поэтому требу- 
ется разработка новых быстрых дискретных ортого- 
нальных преобразований и соответствующих алго- 


ритмов сжатия, пригодных для практической реали- 
зации на современной элементной базе. 

Определим критерии, которым должен соответ- 
ствовать современный алгоритм сжатия изображе- 
ний, пригодный для реализации в аппаратуре. 
Примем размер кадра равным 800x600, количество 
бит на пиксел - 24. 

• Скорость работы. Для современных систем ви- 
деонаблюдения необходима скорость сжатия не 
менее 25 кадров/с. Это означает, что макси- 
мальное время обработки кадра должно быть 
равно 40 мс. 

• Коэффициент сжатия. При высокой скорости 
работы необходимо обеспечить коэффициент 
сжатия изображения не ниже, чем у классиче- 
ского алгоритма ]РЕО, то есть, в 3-25 раз в за- 
висимости от уровня допустимых потерь ин- 
формации. 

• Качество восстановленного изображения. При 
высокой скорости работы необходимо обеспе- 
чить качество изображения, визуально оцени- 
ваемое не ниже, чем у классического алгоритма 
.ІРЕО, в зависимости от уровня допустимых по- 
терь информации. 

• Простота аппаратной реализации. Алгоритмы 
сжатия, работающие в реальном времени, дол- 
жны иметь только две операции - сложение и 
сдвиг, чтобы обеспечить экономию аппаратных 
ресурсов, и, тем самым, их дешевую реализа- 
цию в аппаратном виде, например, с использо- 
ванием недорогих ПЛИС. 

После анализа существующих семейств алго- 
ритмов сжатия изображений с позиции определен- 
ных выше критериев их оценки, авторы пришли к 
следующим выводам: 
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• Алгоритмы сжатия изображений без потерь 
применимы только к узким классам изображе- 
ний (ограниченность цветов, большие области 
одного цвета и т.п.). Неэффективны для сжатия 
полноцветных реалистичных изображений. 

• Вейвлет-алгоритмы (включая стандарт 1РЕС2000) ре- 
курсивны, поэтому для работы требуют много време- 
ни, а также большой емкости запоминающего 
устройства при их аппаратной реализации. Это зна- 
чительно увеличивает стоимость аппаратного кодека. 

• Фрактальное сжатие основано на поиске подоб- 
ных участков изображения, получаемых при по- 
мощи прямого перебора областей изображения 
с применением к этим участкам аффинных пре- 
образований. Даже на современных микропро- 
цессорных платформах программное обеспече- 
ние, реализующее фрактальные алгоритмы сжа- 
тия работает медленно. Сложность аппаратной 
реализации этих алгоритмов также очевидна. 

• Алгоритм сжатия .ІРЕО и его модификации. 
Скорость данных алгоритмов линейно зависи- 
ма от размера изображения, они не требуют от- 
дельного запоминающего устройства при аппа- 
ратной реализации. Также ]РЕО [1] может быть 
усовершенствован, и на основе модифициро- 
ванного алгоритма может быть построен конку- 
рентоспособный кодек. 

Алгоритм .ІРЕО использует дискретно-косинус- 
ное преобразование (Оіхсгеіе-Соаіпе Тгапфтіщ, 
РСТ) [1]. Вычисление ОСТ и обратного преобразо- 
вания ШСТ является одной из наиболее трудоемких 
и длительных процедур, выполняемых при сжатии 
и распаковке кадров. Таким образом, этим процеду- 
ры в наибольшей степени нуждаются в оптимиза- 
ции, позволяющей ускорить работу программного 
и аппаратного обеспечения, производящего сжатие 
и распаковку видеопоследовательностей. 

Строго говоря, дискретно-косинусное преобра- 
зование не является единственно приемлемым ор- 
тогональным преобразованием для сжатия кадров; 
мы можем использовать и любое другое [1]. Рас- 
смотрим подробнее этап вычисления ортогональ- 
ного преобразования. 

Двумерные ортогональные преобразования вы- 
полняются по следующей формуле: 

У = кРХР т . 

Здесь X. 8x8 - матрица входных значений, У: 
8x8 - матрица выходных значений (коэффициен- 
тов), Р: 8x8 - матрица преобразований. А; - скаляр- 
ный масштабирующий коэффициент. Рассмотрим 
теперь матрицу У т : 

У т = (кРХР т ) = кРХ т Р т . 

Введем обозначение Х=РХ Т . Тогда 
У т =к2Р т ; 

У = (к2Р т ) т =кР2 т . 

Таким образом, вычисление матрицы выходных 
значений У может быть сведено к двум последова- 


тельно выполняемым операциям перемножения 
двух матриц, в которых первая матрица одинако- 
вая. Это упрощает программную и аппаратную ре- 
ализацию алгоритма. 

Рассмотрим теперь операцию перемножения 
двух матриц 8x8. Вычисление одного столбца ре- 
зультирующей матрицы в общем случае требует 64 
умножений и 56 сложений. Наши дальнейшие рас- 
суждения будут направлены на то, чтобы умень- 
шить число этих операций для некоторых частных 
случаев ортогональных преобразований (матриц 
Р). Это достигается путем разложения матрицы 
преобразования Дна разреженные множители, т.е. 
представления ее в виде 

Р = ЦхР 2 х...хР т , (1) 

где матрицы Р, содержат максимально возможное 
количество нулей. Такой подход применяется, в 
частности, в 11]. 

Еще одним методом оптимизации двумерного 
ортогонального преобразования является аппрокси- 
мация. Задача аппроксимации сводится к приведе- 
нию коэффициентов матрицы к виду т/2". Это по- 
зволяет избавиться от длительных операций с пла- 
вающей точкой, а также от ресурсоемкой операции 
умножения, заменив ее сложениями и сдвигами. 

Матрица Р для вычисления дискретно-кос- 
инусного преобразования выглядит следующим 
образом [1]: 
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Проведем аппроксимацию преобразования с 
точностью до ІДб- Для этого умножим каждый эл- 
емент матрицы на 16, разделим его на лЙу8=0,3535 
и округлим до четного числа, так, чтобы миними- 
зировать количество получившихся значений. Вы- 
бор числа 16 обусловлен аппаратными ограниче- 
ниями. Дело в том, что для сохранения элемента 
выходной матрицы необходимо 12 бит. Типичная 
разрядность представления таких данных в архи- 
тектуре І386 - 16 бит, то есть, для проведения 
16-битовых операций необходимо умножение каж- 
дого элемента входной матрицы на 2 4 =16. 

После таких преобразований получим матрицу Р\ 
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Естественные науки 


Заметим, что каждая строка матрицы Р' может 
быть представлена в одной из двух следующих форм: 
[а Ъ с (і с? с Ъ а\, 

\а Ъ с й -й -с -Ъ -а]. (2) 

Таким образом, для вычисления произведения 
Р'х Т, где Т=[( 0 ,...,( 1 ] т - вектор-столбец, требуются 
выражения ? 0 -/ 7 ,...,? 4 +? 5 Т 4 — / 5 ; т.е., матрица Р т в раз- 
ложении ( 1 ) имеет вид: 

1 0 0 0 0 0 0 1 

1 0 0 0 0 0 0 -1 

0 1 0 0 0 0 1 0 

0 1 0 0 0 0 -1 0 


0 0 1 0 0 1 0 0 

0 0 1 0 0 -1 0 0 

0 0 0 1 1 0 0 0 

0 0 0 1 -1 0 0 0 


( 3 ) 


Умножение этой матрицы на вектор Т требует I 


Теперь для подсчета Р'х Т требуется 32 умноже- 
ния и 32 сложения, то есть трудоемкость алгоритма 
снижена в 2 раза. Правильность разложения может 
быть проверена перемножением полученных ма- 
триц. Строки 0, 2, 4, 6 матрицы Р ] могут быть пред- 
ставлены в одной из следующих форм: 

[а 0 Ъ 0 Ъ 0 а 0], 

[а 0 Ъ 0 -Ъ 0 -а 0]. 

Используя те же идеи, и опуская промежуточ- 
ные выкладки, получаем: 
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сложений. В итоге получится вектор: 
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Рассмотрим строки, представимые в виде (2). 
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Произведение их на вектор Т имеет вид: 
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и требует 4 умножения и 3 сложения (если уже вы- 
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числен вектор Т). Таким образом, используя (3) и 
(4), получаем разложение матрицы Р'на разрежен- 
ные множители: 
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Таким образом, получение вектора Рх Т потре- 
бует 22 умножения и 16 сложений. Как можно за- 
метить, в качестве множителей для операции умно- 
жения используется всего шесть чисел: 4, 8, 12, 16, 
18 и 22. Приведем операции умножения коэффи- 
циента а на каждый из этих множителей к опера- 
циям сложения и сдвига. 

Исходя из данных табл. 1 , 22 операции умноже- 
ния можно заменить на 14 операций сдвига и 
10 операций сложения. Итого, для получения век- 
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тора Р'х Т необходимо провести 14 операций сдвига 
и 30 операций сложения. То есть, для вычисления 
преобразования для одной матрицы 8x8 необходи- 
мо 2-8Т4 = 224 операции сдвига и 2-8-30=480 опе- 
раций сложения. 


Таблица 1. Замена операций умножения на операции сдвига 
и сложения 


Операция 

Сдвигов 

Сложений 

4а=а<<2=Ь 0 

1 

- 

8а=а<<3=Ь, 

1 

- 

12а=а«3+а<<2=Ь 0 +Ь 1 =Ь 2 

- 

1 

16а=а«4=Ьз 

1 

- 

18а=16а+2а=Ь 3 +а«1=Ь 4 

1 

1 

22а=18а+4а=Ь 4 +Ь 0 =Ь 5 
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Дискретно-косинусное преобразование не яв- 
ляется единственным ортогональным преобразо- 
ванием, которое возможно применить для сжатия 
изображений. Преобразование Уолша [2] характер- 
но тем, что его возможно выполнить при помощи 
одних лишь операций сложения и вычитания. Од- 
нако существует способ дальнейшей оптимизации 
данного преобразования. 

Поставим задачу разложения матрицы преобра- 
зования Уолша Н к на разреженные множители. Для 
этого сначала рассмотрим общий случай - блоч- 


ную матрицу 


Т = 


О О 

О -о. 


, где 0 - произвольная 


матрица размерности тхт. Разложение матрицы Т 
на разреженные множители имеет вид: 
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где Е т - единичная матрица размерности тхт , 
О т - нулевая матрица размерности тхт. 

Действительно, пусть Х=(х ь ...,х 1т ) Т - вектор раз- 
мерности 2/я, Х'=(х ь ...,х т ) т , Х"=(х т+ь -,х 2т ) Т , тогда: 
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С другой стороны, 
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(5) 


Таким образом, разложение верно. Вернемся 
теперь к матрице Н к . Для разложения ее на разре- 
женные множители, уточним ее структуру. Введем 
матрицы: 



"1 Г 


~А 

А 

А = 

1 -1 

, В = 

А 

-А 


С использованием этих обозначений матрицу 
Н к можно записать: 


Я 


И 


в в 
в -в 


Используя (5) и (6) получаем: 
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Таким образом, матрица Н к может быть пред- 
ставлена как произведение трех разреженных мно- 
жителей: 


Я 


И 


1 0 0 0 1 0 0 0 

0 1 0 0 0 1 0 0 

0 0 1 0 0 0 1 0 

0 0 0 1 0 0 0 1 

1 0 0 0 -1 о о о 

0 1 0 0 0 -1 о о 

0 0 1 0 0 0 -1 о 

0 0 0 1 0 0 0 -1 


1 0 1 

0 1 о 

1 0 -1 

0 1 о 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

"1 1 о 

1 -1 о 

0 0 1 

0 0 1 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


0 0 0 0 0 
1 0 0 0 0 
0 0 0 0 0 
-10 0 о о 
0 10 10 
0 0 10 1 
0 10-1 о 
0 0 1 0-1 

0 0 0 0 о 

0 0 0 0 0 

1 0 0 0 0 

-1 0 0 0 о 

0 110 0 
0 1-10 о 
0 0 0 1 1 
0 0 0 1 -1 


Для умножения матрицы преобразования на 
один вектор-столбец исходных данных необходимо 
всего 24 сложения. Это означает, что для преобразо- 
вания всей матрицы 8x8 необходимо 2-8-24=384 
сложения, т.к. всего вектор-столбцов в матрице ис- 
ходных данных восемь, и для вычисления двумерно- 
го преобразования матричные умножения на разре- 
женные множители необходимо производить дваж- 
ды - для вектор-строк и для вектор-столбцов. 

Тестирование разработанного алгоритма прово- 
дилось на компьютере Репйшп IV - 1,7 ОН 2 / 256 МЬ 
ООП ВАМ / 120 СЬ ИБО, сііірчеі Іпіеі 1845 О под упра- 
влением ОС МІСГО 80 Й ЧѴіпсіоѵл ХР (§егѵісе Раск 2). Ис- 
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Таблица 2. Результаты тестирования алгоритмов, построенных на основе быстрого и классического преобразования Уолша 


Имя файла 

Размер 
файла ис- 
ходный, 
байт 

Размер 

сжатого 

файла, 

байт 

Кс* 

Быстрое двумерное преобразование Уолша 

Классическое преобразование Уолша | 

Тех, МС 

Трасп # МС 

Скорость 

сжатия, 

кб/с 

Скорость 

расп-ки, 

кб/с 

Тс*, мс 

7расп, МС 

Скорость 

сжатия, 

кб/с 

Скорость 

расп-ки, 

кб/с 

сІедд.Ьтр 

2148978 

756439 

2,84 

1390 

781 

1546,03 

2751,57 

1610 

1344 

1334,77 

1598,94 

Ігутіге.Ьтр 

3706188 

1192680 

3,11 

1859 

1406 

1993,65 

2635,98 

2828 

2328 

1310,53 

1592,01 

Іепа.Ьтр 

786450 

222854 

3,53 

391 

313 

2011,38 

2512,62 

625 

531 

1258,32 

1481,07 

топагсЬ.Ьтр 

1179666 

288107 

4,09 

672 

656 

1755,46 

1798,27 

859 

735 

1373,30 

1604,99 

реррегз.Ьтр 

786450 

196604 

4,00 

500 

500 

1572,90 

1572,90 

547 

640 

1437,75 

1228,83 

заіІ.Ьтр 

1179666 

401513 

2,94 

1094 

531 

1078,31 

2221,59 

1313 

812 

898,45 

1452,79 

зеггапо.Ьтр 

1498296 

415437 

3,61 

953 

578 

1572,19 

2592,21 

1047 

1094 

1431,04 

1369,56 

ІиІірз.Ьтр 

1179666 

411649 

2,87 

875 

500 

1348,19 

2359,33 

969 

797 

1217,41 

1480,13 

Итого: 

12465360 

3885283 

3,21 

7734 

5265 

1611,76 

2367,59 

9798 

8281 

1272,24 

1505,30 


Примечание: К а - коэффициент сжатия, Т^, Г расп - время сжатия и распаковки 


Таблица 3. Результаты тестирования алгоритмов, построенных на основе аппроксимированного и классического преобразова- 
ния ОСТ 


Имя файла 

Размер фай- 
ла ИСХОД- 
НЫЙ, байт 

Аппроксимированное ОСТ 

Классическое ЭСТ 

Кех 

Тс ж, мс 

Траст. МС 

Скорость 
сжатия, кб/с 

Скорость 
расп-ки, кб/с 

Ксх 

Тс ж, мс 

Траст/ МС 

Скорость 
сжатия, кб/с 

Скорость 
расп-ки, кб/с 

сІедд.Ьтр 

2148978 

2,64 

1441 

954 

1491,20 

2251,55 

2,16 

2594 

1718 

828,44 

1250,86 

Ігутіге.Ьтр 

3706188 

3,07 

2422 

1675 

1530,08 

2212,65 

2,51 

4360 

3015 

850,04 

1229,25 

Іепа.Ьтр 

786450 

3,37 

495 

373 

1588,79 

2106,56 

2,76 

891 

672 

882,66 

1170,31 

топагсЬ.Ьтр 

1179666 

4,41 

712 

556 

1657,61 

2123,40 

3,62 

1281 

1000 

920,89 

1179,67 

реррегз.Ьтр 

786450 

4,12 

503 

391 

1562,48 

2013,67 

3,38 

906 

703 

868,05 

1118,71 

заіІ.Ьтр 

1179666 

3,06 

929 

556 

1269,98 

2123,40 

2,51 

1672 

1000 

705,54 

1179,67 

зеггапо.Ьтр 

1498296 

3,49 

894 

677 

1675,11 

2214,23 

2,86 

1610 

1218 

930,62 

1230,13 

ГиІірз.Ьтр 

1179666 

2,86 

990 

521 

1191,58 

2266,17 

2,35 

1782 

937 

661,99 

1258,98 

Итого: 

12465360 

3,16 

8387 

5702 

1486,33 

2186,27 

2,59 

15096 

10263 

825,74 

1214,59 


пользовался международный стандартный тестовый 
набор Шііегіоо Вщіопе ( СаІ§агу Согрш ) [3] и специаль- 
но разработанный программно-экспериментальный 
комплекс [4]. Результаты тестирования алгоритмов 
сжатия на основе рассмотренных ортогональных пре- 
образований приведены в табл. 2 и 3. 

Разработанный метод вычисления преобразо- 
вания Уолша при том же качестве восстановленно- 
го изображения превосходит классический метод 
вычисления дискретно-косинусного преобразова- 
ния по скорости сжатия и распаковки в 1,9 раза. 

Разработанный аппроксимированный метод вы- 
числения дискретно-косинусного преобразования 
при том же качестве восстановленного изображения 
превосходит классический метод вычисления дис- 
кретно-косинусного преобразования по скорости 
сжатия и распаковки в 1,8 раза. Сравнительная тру- 
доемкость преобразований приведена в табл. 4. 

Коэффициент сжатия для обоих преобразований 
возрастает примерно на 22 % по сравнению с класси- 
ческим дискретно-косинусным преобразованием. 

Итак, авторами разработаны и протестированы 
упрощенный метод вычисления дискретно-кос- 
инусного преобразования и быстрое двумерное 
преобразование Уолша. Оба они превосходят клас- 
сический метод ОСТ более чем в два раза и по сжа- 
тию, и по распаковке, при том же, или более высо- 
ком коэффициенте сжатия и высоком качестве 
восстановленного изображения. 


Таблица 4. Трудоемкость преобразований ( число операций 
для матрицы 8x8) 


Преобразование 

Сложений 

Сдвигов 

Умножений 

Классическое ОСТ 

896 

- 

1024 

Аппроксимированное ОСТ 

480 

224 

- 

Классическое преобразова- 
ние Уолша 

896 

- 

- 

Быстрое двумерное преобра- 
зование Уолша 

384 

- 

- 


Для оценки качества изображения используется 
такая величина, как среднеквадратичное отклоне- 
ние восстановленного изображения от исходного. 


Зависимость среднеквадратичного отклонения от 
коэффициента сжатия 



— - — Среднеквадратичное отклонение 4РЕО 
^^^^“Среднеквадратичное отклонение ѴѴНТ 


Рис. 1. Качество восстановленного изображения 

На рис. 1 приведен график зависимости сред- 
неквадратичного отклонения от коэффициента 
сжатия. На рис. 2-5 можно увидеть исходное изо- 
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Рис 2. Исходное 
изображение 


Рис.З. Изображение, 
восстановленное 
после сжатия с ис- 
пользованием 
классического 
ОСТ при К сж =2,76 


Рис. 4. Изображение, 
восстановленное 
после сжатия с 
использованием 
аппроксимиро- 
ванного ОСТ при 
Ксж=3,37 


Рис. 5. Изображение, 
восстановленное 
после сжатия с 
использованием 
быстрого преоб- 
разования Уолша 
при /Сок =3,53 


бражение и результаты работы алгоритмов сжатия 
и распаковки с использованием рассмотренных 
выше преобразований. 

Упрощенный метод вычисления дискретно- 
косинусного преобразования совместим с любым 
другим методом сжатия изображений, имеющим в 
своей основе дискретно-косинусное преобразова- 
ние (]РЕО, МРЕС ), поэтому его очень удобно ис- 
пользовать для быстрой распаковки уже имеющих- 
ся сжатых изображений или видеопоследователь- 
ностей. Поскольку данное преобразование факти- 
чески является разновидностью ЪСТ , средний ко- 
эффициент сжатия алгоритма на базе этого преоб- 
разования сравним с коэффициентом сжатия алго- 
ритма .ІРЕСі. Рекомендуемые сферы применения: 


программные и недорогие аппаратные кодеки (на- 
пример, в цифровых фотоаппаратах) для записи, 
хранения и просмотра видеопоследовательностей в 
широкораспространенных форматах. 

Быстрое двумерное преобразование Уолша име- 
ет низкие требования к аппаратуре и высокую ско- 
рость работы при достаточно высоком качестве 
восстановленного изображения, хотя и несколько 
уступающем алгоритму ,/РЕО. Это делает данное 
преобразование эффективным применительно к 
программным и аппаратным системам видеона- 
блюдения, видеосвязи и видеоконференций, си- 
стемам терминального доступа, а также ко всем 
приложениям, требующим быстрой передачи ви- 
деоданных. 
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